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Abstract 

We study the individual contribution to secondary lepton production in hadronic interactions of 
cosmic rays (CRs) including resonances and heavier secondaries. For this purpose we use the same 
methodology discussed earlier [Ij], namely the Monte Carlo particle collision code DPMJET3.04 to 
determine the multiplicity spectra of various secondary particles with leptons as the final decay 
states, that result from inelastic collisions of cosmic-ray protons and Helium nuclei with the in- 
terstellar medium of standard composition. By combining the simulation results with parametric 
models for secondary particle (with resonances included) for incident cosmic-ray energies below a 
few GeV, where DPMJET appears unreliable, we thus derive production matrices for all stable 
secondary particles in cosmic-ray interactions with energies up to about 10 PeV. 

We apply the production matrices to calculate the radio synchrotron radiation of secondary 
electrons in a young shell-type SNR, RX J1713. 7-3946, which is a measure of the age, the spectral 
index of hadronic cosmic rays, and most importantly the magnetic field strength. We find that the 
multi-mG fields recently invoked to explain the X-ray flux variations are unlikely to extend over a 
large fraction of the radio-emitting region, otherwise the spectrum of hadronic cosmic rays in the 
energy window 0.1-100 GeV must be unusually hard. 

We also use the production matrices to calculate the muon event rate in an IceCube-like detector 
that are induced by muon neutrinos from high-energy 7-ray sources such as RX J1713. 7-3946, Vela 
Jr. and MGRO J2019+37. At muon energies of a few TeV, or in other word, about 10 TeV neutrino 
energy, an accumulation of data over about flve to ten years would allow testing the hadronic origin 
of TeV 7-rays. 

PACS numbers: 13.85.Tp, 95.85.Ry, 98.70.Sa, 98.38.Mz 

Keywords: cosmic rays, cosmic-ray interactions, neutrino and lepton astronomy, supernova remnants 



I. INTRODUCTION 

It is known that hadronic cosmic rays interacting with the interstellar medium (ISM) can 
produce 7-rays and leptons. The principal production mechanisms for 7-rays in high-energy 
astrophysics are inelastic CR+ISM interactions with subsequent decays of the secondaries 
(mostly the neutral vr" and 77) into 7-rays. The decays of charged vr^ mesons and other 
secondaries will result in the production of high-energy neutrinos and secondary leptons, 
electrons and positrons. In GeV-scale Galactic cosmic rays, the secondary electrons and 
positrons comprise around 20% of the total electron flux and significantly contribute to the 
non-thermal electromagnetic radiation of 7-ray sources. The flux of cosmic-ray neutrinos 
is determined by the cosmic-ray abundance and the production cross sections of the par- 
ent particles (vr^, i^^, K^, i^£, S^, S°, A, A) in the collisions of cosmic-rays with the 
interstellar gas nuclei. 

The purpose of this work is to carefully calculate the lepton production in cosmic-ray 
interactions. For this purpose we use the high-energy physics event generator DPMJET-III 
|2j to simulate all secondary productions in both p-generated and He-generated interactions, 
similar to earlier work on hadronic 7-ray production pj. Our study includes direct lepton 
production and the decays of all relevant secondary particles with leptons as the final decay 
products. For the cosmic rays, protons and helium nuclei are taken into account. We 
assume the composition of the ISM as 90% protons, 10% helium nuclei, 0.02% carbon, and 
0.04% oxygen. At energies below 20 GeV, where DPMJET appears unreliable, we combine 
the simulation results with parametric models for vr^ and 7r° production |3|, |j] that include 
the production of the resonances A(1232) and A(1600) and their subsequent decays, thus 
deriving a lepton production matrix for cosmic rays with energies up to about 10 PeV that 
can be easily used to interpret the spectra of cosmic lepton sources. 

We use our production matrices to calculate the CR-induced electron/positron production 
in young shell-type supernova remnants (SNR), for which the neutrino fluxes are also cal- 



culated. Observations of non-therma 



RX J1713. 7-3946 



g, IC 443 



X-ray synchrotron radiation from the SNRs SN 1006 



Gas A |9|, and RGW 86 [lO| support the hypothesis 
that the Galactic cosmic-ray electrons may be accelerated predominantly in SNR, although 
other acceleration sites of high-energy electrons might exist in the Galaxy. The existence of 



electrons in SNRs with an energy up to about 100 TeV implies TeV-scale 7-ray emission 



Q 



besides that possibly produced by hadronic cosmic rays, and TeV-scale gamma-rays have 



indeed been detected from the SNRs Vela Junior [12] and RX J1713. 7-3946 [JJ]. Although 



the hadronic interpretation is often favored to explain the 7-ray emission from SNRs such as 
RX J1713. 7-3946 [l3|, there is still no direct observational evidence for nucleon acceleration 
in SNRs pj, uJ^. Currently, both hadronic and leptonic models can describe the observed 



TeV 7-rays in different astrophysical objects, but have problems of similar magnitude 15 1. 



RX Jl 713. 7-3946 is a unique SNR in the sense that its X-ray emission is strongly domi- 
nated by a non-thermal component, which is p resumed to be synchrotron radiation of ultra- 
relativistic electrons (see references Il6l.ll7l.ll8ll and references therein for details). The recen; 

n r 

broadband X-ray spectroscopy performed with Chandra [I8j and the Suzaku experiment [19 1 
provides evidence of very effective acceleration of particles in the shell of RX J1713. 7-3946 
and also reveals that the X-ray emission from two compact regions is variable in flux, which 
the authors interpret as indication of multi-mG magnetic fields. If this conclusion is correct 
and applies to the entire remnant, then the observed TeV-band 7-ray emission must be of 
hadronic origin. Here we show that the synchrotron flux from the secondary electrons in 
RX J1713. 7-3946 would in fact exceed the observed radio flux from this object, if multi-mG 
magnetic fields would permeate the remnant, so a magnetic field of such magnitude will 
likely exist at most in a small fraction of emission region. 

In a hadron accelerator, TeV neutrinos should be produced in roughly the same number 
as TeV 7-rays. The Universe is generally more transparent for high-energy neutrinos than 
for TeV and PeV 7-rays, that can be absorbed by pair production with either the microwave 
background (CMB) or the infrared/optical photon background (IRB). Therefore the detec- 
tion of astrophysical neutrinos can provide invaluable complementary information on the 
existence of hadronic processes at astrophysical objects. Here we calculate the observable 
TeV-scale neutrino flux from RX J1713. 7-3946 and other sources. This is not the first at- 
tempt to do such work. In fact, several calculations of the expected 7-ray and neutrino 
spectra were made for parametrized hadronic particles in the energy range of neutrino de- 



tectors 
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26(1 . In contrast to the earlier works, we include all relevant 
channels of neutrino production and also we do not rely on assuming a simple power-law 
proton spectrum or a power-law with an exponential cutoff, because our neutrino production 
matrix can be applied to energetic particles with arbitrary spectrum. 



II. INELASTIC COSMIC-RAY INTERACTIONS 

A. Set-up of the Monte-Carlo event generator DPMJET-III 

This work uses the same method to determine the 7-ray production matrix for the cosmic- 
.ay had.onic inte.act.ons as published earfie. fo. .-.ay p.oduct.on fl, na„.e.y apply.ng the 
Monte-Carlo event generator, DPMJET-III [2|], to simulate the secondary production in 
cosmic-ray interactions. The readers are referred to reference [l| for the details of this 
technique. 

B. Decay channels of secondary particles to leptons 

All secondary particles produced in simulated hadronic interactions are recorded while 
running the event generator DPMJET-III. The following decay modes show all the decay 
processes of secondary products with leptons as the final decay particles, that are taken 
into account. Stable leptons considered in this work are e^, z/g, t'e, t'/x, and z/^. Note that 
the lifetime of neutrons is about 886 s, much shorter than the propagation time scale of 
cosmic-ray particles in the Galaxy. Therefore, the neutrons are treated as having decayed 
entirely in this work. The decay channels considered in this work are as follows. 



baryonic decays: 


n - 


^p + e + z/e. 


n - 


■^ p + e+ + Ue, 


A- 


ip+--, 




[n + vrO, 


A- 
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p + vr+, 
-> < 




ra + 7r°. 


EO 


^ A + 7, 


S+ 


^ f p + VfO, 




1 n + 7r+. 


s- 


^ n + TT". 


mesonic 


decays: 




leptonic decays: 



/i"^ ^ e+ + Z/e + z/^, 
jj,~ ^ e" + Ue + y^, 

e+ + Z/e + z/^, 

e" + z/g + z/^, 

In the DPMJET simulation setup, the decay of r] mesons produced in the hadronic 
interactions has aheady been treated in the simulation. In our study on cosmic-ray induced 
7-ray production |l|], we have verified that the rf mesons are properly accounted for by 
performing test runs under different PYTHIA parameters, for which the rj decays can be 
avoided in the simulation. We have also performed another independent verification using 
the very small contribution to e~ /e^ and /i^/;U+ pairs from r] decays, which are found 
consistent within the statistical uncertainty with all 7] mesons having decayed. 



The TT^ decays and subsequent /i^ decays via the it — ^ — e decay chain are dominant 
in cosmic-ray interactions because pions carry the highest multiphcity among all the sec- 
ondaries. We ignore neutrino production by charmed particles such as D and D, because 
they contribute only at E^, > 100 TeV 27|], where they are ignorable compared with neu- 
trino production by pions. In the case of decays of moving vr^'s, the vr energy is roughly 
equally divided among the four decay products, and thus the flux ratios are approximately 
(z/g + i'e)/{i'^ + Vy) = 1/2 and v^^/v^ = 1, no matter what the vr spectrum is. In reality, 
however, these neutrino ratios may vary on account of the decays of heavier baryons and 
mesons. An excess of /i"*" multiplicity arises because more 7r"^'s than 7r~'s are produced, and 
so the ratios I'e/'^e and v^/v^ are somewhat increased. 

For the two-body decay processes, the decay spectra are evaluated by particle kinematics; 
for the electrons/positrons from muon decays, the decay spectra are calculated by Lorentz 



transformation of the particle distribution in the center-of-mass system of the muon J28|, 



which also includes the effect of the /z polarization [29|, |30|; for kaon decays and heavy 
nucleon decays, the Dalitz-plot distribution is applied. The values published in Particle 
Data Group are employed for the fraction of each individual decay process. 

C. Resonance contribution at low energies 

In our earlier work on 7-ray production \v\ we already found that DPMJET is unreliable 

)elow a few GeV collision energy. At these energies we therefore use a parametric model 

^ to calculate the resultant decay spectra of 7, e^, v^ {^e) and z/^ (z/^). For the decay 

of the resonances A(1232) and A(1600) we assume the pion momentum to be isotropically 

distributed in the center-of-mass system and no angular correlation between the two pions 

in A(1600) decay. 

Note that the parametric model for pion and resonance production at low incident energies 
is for pp interactions only. We then use DPMJET to calculate the energy-dependent weight 
factors, which allow us to parametrically account for p+ISM and He+ISM collisions in the 
parametrization approach, even though it is derived for pp collisions only. The weight factors 
carry no strong dependence on the energy of the projectile particle. 



III. THE PRODUCTION MATRIX OF SECONDARY LEPTONS GENERATED 
BY COSMIC RAYS 

The differential production rate of a final secondary particle is given by 

U/Ti f ( dji \ 

'^-^f^) = MriETW = "'-' A,. ''^™ '^-<^-«) '^-- ["Ie) ('' 

with Ncr{Ecr) = dE^^'^dv (GfeV cm )^^ as the differential density of CR particles (p or a). 
The differential cross section of a secondary particle produced in (p,Q;)-ISM collisions is 

da . , dn , . 

— {Ecn,E) = a,rodj^ (2) 

where aprod is the inelastic production cross section, i.e., the sum of diffraction, non- 
diffraction and resonance components whose values are calculated by DPMJET and the 
parametric model; ^ is the multiplicity spectrum of the secondary particle in question. 

We follow the decay of unstable secondary particles to stable particles, and thus obtain 
the final spectrum of the stable secondaries 7-rays, e^, Ue, Ue, u^ and z/^ as 

Q2ndiE.d = V njsM [ dEcR Ncr{Ecr) cf3cR cj{Ecr) ^(^i, Ecr) (3) 

fc J Ecr "-^j 

where -j^ is the multiplicity spectrum of stable particle species i resulting from production 
channel fc, either an unstable secondary or direct production. For binned particle spectra 
the production integral, Eq. ([3]), can be re-written as 

Q2UE,) = njsM J2 ^^^ ^cr{E,) ci5, aiE,) ^ ^(E„ E,) (4) 



dEi 

j k 



nisM 



J2^E,NcR{E,)cf3,ajM,, (5) 



J 



for the secondary particle of interest. In Eqs. (jlj) and (^, the energies for cosmic-ray 
particles {p or a) and the secondary particles are parametrized as [jj: 

Et = 1.24 ■ (1 + 0.05)^' GeV/n, for cosmic-rays, (6) 

Ek = 0.01 ■ (1.121376)^"°-^ GeV, for secondary particles. (7) 

nisM — 1-11 n-H in Eq. (^ accounts for the specific element abundance in the ISM assumed 
90% H, 10% He, 0.04% O and 0.02% C here. Please note that the ratio of the 7-ray and 
neutrino yields is indepedendent of ujsm- The main product of our work is the production 
matrix Mjj that describes the production spectrum of the stable secondary particles, one 
for each, for arbitrary cosmic-ray spectra, separately for protons and Helium nuclei. 

8 



IV. APPLICATIONS: SNR 

With the production matrix estabhshed for each stable particle, we are now in the po- 
sition to calculate the production spectra of secondary leptons and neutrinos in a variety 
of sources. As an example we here consider young shell-type SNRs. Being associated with 
dense gas along the line-of sight, the shell-type SNR RX J1713. 7-3946 could provide signifi- 
cant information on the acceleration of hadronic cosmic-rays. In addition, RX J1713. 7-3946 
is also a unique supernova remnant in the sense that its X-ray emission is strongly dominated 
by a non-thermal component, which is presumed to be synchrotron radiation of ultrarela- 
tivistic electrons. These electrons also produce TeV-band gamma-rays, but would require a 
relatively small magnetic field to account for the TeV-band 7-ray spectrum observed with 



HESS 131]. The recently observed rapid variability of X-rays of RX J1713. 7-3946 has been 



interpreted as evidence of a multi-mG magnetic field 18| , in which case the leptonic scenario 
would be untenable. 

We can test this interpretation by calculating the synchrotron radiation of secondary 
electrons in this remnant. In an earlier paper we have determined the cosmic-ray spectrum in 
RX Jl 713. 7-3946 under the assumption that the observed TeV-band spectrum is hadronic in 
origin [l| . Using our newly derived production matrix for secondary electrons and positrons, 
we can determine the differential production rate of those particles which, together with a 
time-dependent cosmic-ray continuity equation, allows us to calculate the radio synchrotron 
radiation of the secondary electrons with mainly the magnetic field strength as a fully free 
parameter. 

Our fit of the TeV 7-ray spectrum constrains only the strongly curved cosmic-ray spec- 
trum in the 1-100 TeV energy range. However, synchrotron emission in the GHz range is 
produced by GeV-scale electrons and we do not know how the spectrum continues down to 
GeV energies. So for simplicity, we match a single-power-law spectrum to the multi-TeV 
spectrum of cosmic-ray nucleons that we derive from fitting the HESS data, 

N{E) = No i^—j exp [-^^J , E>E, (8) 

N{E) = Mo (■^] , E<E, (9) 

where the parameters in equationJSj) are derived by fitting the observed TeV-band 7-ray 
spectrum of RX J1713. 7-3946 l|, |3l|- In this parametrization, Eq = 15 TeV is a normaliza- 



tion chosen to render uncorrelated the variations in the power-law index, s, and the spectral 
curvature, a. -Emax is the cut-off energy, and the best-fitting parameter values are s = 2.13, 
a = —0.25, and -Emax ^ 200 TeV. iVn is the overall normalization, determined by the integral 



spectrum of 7-rays above 1 TeV 13|. If a numerical value of a is set, the normalization 
Mq and the merging energy Ec are fully determined by the two continuity conditions that 
link the two spectral forms ([8]) and iQ. Models of particle acceleration in cosmic-ray mod- 



ified shocks suggest that the spectral index a should be slightly smaller than 2 [32|, 
Note that a denotes the average particle spectral index between about 1 GeV and 10 TeV. 
The combined cosmic-ray spectrum together with our production matrices for electrons and 
positrons then yields the differential production rate of secondary leptons that is implied by 
the 7-ray spectrum observed from RX J1713. 7-3946. 

A. Electrons/Positrons 

Secondary electrons have been accumulated in the remnant since it commenced hadron 
acceleration. To determine the synchrotron radiation of the secondary electrons, we must 
obtain the electron spectrum in the SNR by solving a simplified time-dependent electron 
transport equation: 

^ + 4(Ki^)A^) = Q(^,2^) (10) 

where T is the time after supernova explosion and the electron energy loss rate is that for 
synchrotron radiation 

6(E) = E = -^aThct/B/32 7' = -^^ThC^/3'7' (H) 

with Ub as the energy density of the magnetic field and a^h as the Thomson cross section. 
The likely distance to RX J1713. 7-3946 is 1 kpc and the implied age would be about 



1,600 years p. 



18 



34i . An alternative kinematic solution places the remnant at about 6 kpc 

n 

in distance with implied age of 10,000 years j8l]. Whatever the acceleration history of the 
remnant is, cosmic-ray particles which have been accelerated early in its evolution will have 
lost some of their energy by adiabatic expansion, with most of it happening very early in 
the evolution. To be conservative we therefore consider only half of the remnant age, during 
which we assume the differential production rate of secondary electrons, q{E), to be constant 

10 



and given by our production matrix and the cosmic-ray spectrum as in Eqs. ([8]) and iQ 

Q{E,T) = q{E)Q{T-ty2) . (12) 

Then the spectrum of secondary electrons or positrons is obtained by integrating the Green 
function and the source term over time and energy: 

N{E, T) = J J dEo dTo q{E,) Q{T, - ti/2) ^^^ 7Jl7 ^^ ® (^0 " ^) (l^) 

where t = J^ -4^ is the electron cooling time. For high energies r < ti/2 and the spectrum 
is loss-limited, i.e. essentially steepened by 1 in the spectral index, whereas at low energies 
the production spectrum is preserved. Note that the synchrotron flux calculated using the 
electron spectrum ( iT3ll is absolutely normalized by the TeV 7-ray flux. 

Figure [1] shows an example of the synchrotron spectrum from RX J1713. 7-3946 that is 
contributed by secondary electrons assuming a distance of 1 kpc, the age as 1,600 years, 
and the power-law index a = 1.8 for the multi-GeV cosmic-ray nucleon spectrum. The 
synchrotron spectra for three different values of the magnetic fields, B = 500, 2000, and 6000 



/iC, are shown in comparison with the observed X-ray [16] and radio 13|, ISSj data. At higher 
frequencies the synchrotron power is always about half the TeV-band 7-ray power on account 
of the equal likelihood of neutral and charged pions being produced in hadronic collisions. 
The non-thermal X-ray emission is therefore always contributed by primary electrons. 

To be noted from the figure is that for a field strength of 2 mG the radio synchrotron flux 
is similar to the observed flux, so the primary electrons must be very few because they would 
contribute only the remaining fraction of the observed radio flux, thus further lowering the 
e/p ratio in accelerated particles [15]. This limit is much more severe, if the distance to 
RX J1713. 7-3946 is more than 1 kpc, and the age correspondingly more than 1,600 years. 
It is also more severe, if the effective GeV-to-TeV spectral index of cosmic-ray nucleons in 
the source is softer, a > 1.8, and it is relaxed in the opposite case. One should note that 
in mG-fields GHz-band synchrotron emission is emitted by electrons with about 300 MeV 
kinetic energy, i.e. below the relativistic transition. The effective spectral index must be 
measured by comparing the 1 GeV flux and the ~ 10 TeV flux of cosmic-ray nucleons, which 
is considerably softer than the local spectral index near 1 TeV. Measurements of the radio 
synchrotron of a large number of shell-type SNRs indicate that in the GeV band, the spectral 



index is locally around 2 |36|, whereas at 1-10 TeV the local spectral index should be closer 

11 



to 1.5 [32], so on average we may expect a ^ 1.8. Consequently, a multi-mG magnetic 
field is unlikely to exist in a large fraction of the synchrotron emission region or a must 
be significantly smaller than 1.8 (i.e., the spectrum is harder). This result is conservative, 
because we are ignoring secondary electrons that were produced in the first half of the time 
since the SN explosion. 

The recent observation of rapid variability of X-rays of RX J1713. 7-3946 [18] has been 
interpreted as evidence for the notion that the magnetic field of the hot spots (two compact 
regions) could be as large as mG level. These compact regions are in only a few arcseconds 
and likely present transient acceleration sites of high-energy electrons. Our calculations 
shows that, if such strong magnetic field exists in SNR shell, they are unlikely to fill the 
entire remnants. 

Table 1 summarizes the limits on the averaged magnetic-field strength for various values 
of the effective GeV-to-TeV spectral index of cosmic-ray nucleons and different values for the 
age of RX J1713. 7-3946. The contribution of synchrotron radiation by secondary electrons 
is compared with the integrated fiux density 20 Jy at 1.4 GHz, which is the result of 35] 
for the NW region but scaled to the entire remnant with a factor 3 based on the angular 
distribution of TeV gamma-ray emission 13| . This scaling factor carries a sizable uncertainty 
not less than 30%, largely owing to systematic uncertainty in the background subtraction in 



the radio data. Some similar works 
remnant. 



13 



35l | have used a scaling factor of 2 for the entire 



B. Neutrinos 

High-energy cosmic neutrinos are inevitably produced in parallel with hadronic 7-rays. 
Therefore, TeV-band 7-ray sources represent the prime targets to search for cosmologi- 
cal neutrinos. Several calculations on the expected 7-ray and neutrino spectra were made 



for parametrized spectra of hadronic 



E 



1000) TeV 
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partic 
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es in the energy range of neutrino detectors, 



261], usually by assuming a simple power-law 



proton spectrum or a power-law with an exponential cutoff. In contrast, our neutrino pro- 
duction matrix can be applied to cosmic-ray nucleon spectra of any form. Also, the particle 
yields in the parameterizations can be underestimated because they are often based only 
on the pion channels. Our production matrix accounts for all relevant decay processes in 
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cosmic-ray hadronic interactions. As shown in Figure [2], we find that using the same cosmic- 
ray spectrum, our production matrix gives about 15-20% more 7-rays at all energies. Our 
production matrix also yields an approximately 30% more muon neutrinos than the para- 
metric model 25| after full mixing. The excess of 7-rays and neutrinos in this figure agrees 
with the conclusion of our earlier study 11], which has considered the full picture of 7-ray 
production in hadronic interactions and indicated about 20% more hadronic 7-rays relative 
to 7-rays from vr*^ decays alone. Also note that ?7-decays contribute about 8-10% and kaons 
also have roughly the same contribution [l| . Only a few percent of 7-rays are contributed 
by other decay processes (together with the direct 7-ray production). For a given 7-ray or 
neutrino flux, our production matrix therefore implies a lower flux of cosmic-ray nucleons 
than estimated using the parametric models 2J, |25 |. 

In high-energy hadronic interactions, the resulting initial v flavor ratio from vr decays is 
Vf^ : Vy^ : Vt = \ : 2 : ^. However, the neutrino oscillations transform this ratio to 1:1:1, 
i.e., full mixing, because the SNR size of ~ 10 pc implies a substantial range of neutrino 
path lengths, wider than the neutrino oscillation wavelength Lose = (47r£'/Am^)^c ~ 
0.8 (i5/TeV)(Am^/10~^°eV )^^ pc [37]. Therefore, it is reasonable to expect the v^ spectrum 
arriving at Earth to be about 1/3 of the total dN^/dE. The resulting v^ spectrum is o nly 



weakly dependent of the relative fraction of z/g's and i/^'s at the source. Other studies |21l. |22] 
have calculated high-energy neutrinos from galactic sources considering the actual neutrino 
mixing angles for point sources. For a comparison of the neutrino source functions by our 
production matrix and the parametric model 221], we therefore must consider the v^ spectra 
from the TeV 7-ray source RX J1713. 7-3946 before neutrino oscillations. In Figure Ej the 
v^j_ spectrum is calculated using the production matrix and a generating cosmic-ray particle 
spectrum 1| that is determined as best-fltting the observed HESS 7-ray data of RX J1713.7- 
3946; the v^ spectrum by the parametric model [22| is calculated based on the same 7-ray 
spectrum. As seen in this flgure, the parametric model gives overall lower raw z/^ yield from 
RX J1713.7-3946. 

Practically, neutrinos and anti-neutrinos are not distinguishable in neutrino telescopes. 
Therefore, full neutrino mixing is usually assumed for the v spectrum. Figure [3] shows the z/^- 



induced muons which may be detected in an IceCube-like experiment [26 
sources such as RX J1713. 7-3649 [jj, MGRO J2019+37 ^ and Vela Jr. 



from TeV 7-rays 



381 ] , together with 



the background from atmospheric neutrinos. The z/^-induced muons for each source example 
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and the atmospheric background are calculated by a parametric model 39|, including the 
neutrino attenuation due to scatterings within the Earth. 

In calculating the i/^-induced muon event rate, the generating cosmic-ray particle spec- 
trum for each source was obtained by fitting the observed 7-ray spectra using our gamma-ray 
production matrix l|], from which the i^-spectra are determined using the neutrino produc- 
tion matrix ([5]) and assuming full mixing. Note that a TeV 7-ray spectrum from MGRO 
J2019+37 is not available at this time. Therefore a fit was performed using EGRET data 
for 3EG J2021+3716 26|. A more reliable calculation will depend on a measurement of the 
TeV 7-ray spectrum of this source. 

The charged-current and neutral-current cross sections, the inelasticity involved in vN — >■ 
IM^ + X interactions, and the range of the muons in the detector are based on the theoretical 



results of 



39 



40 



4l|. The //-induced muon event rate for the atmospheric background is 



then calculated using the parametrized atmospheric neutrino background as derived in 42| , 
both vertically and horizontally, and assuming an angular resolution element of 3 deg^ (1 
deg uncertainty radius), which roughly corresponds to the angular resolution of IceCube. As 
seen in Figure [3l the calculations indicate a good detectability of high-energy neutrinos from 
the high-energy 7-ray sources at a few TeV in muon energy, or about 10 TeV neutrino energy. 
An accumulation of data over about five to ten years would allow testing the hadronic origin 
of TeV 7-rays. 

The lepton and neutrino production matrices are applicable for arbitrary cosmic-ray 
hadron spectra. The main systematic uncertainty is that of the simulation tool DPMJET 



itself. DPMJET has been extensively tested against data 
43|, for the simulation of cosmic-ray air shower problems 



or accelerator experiments 



m 



4J, |45|] , and for the atmospheric 



7-ray and neutrino production 46|. Although the uncertainty for hadronic productions 



simulated by DPMJET differs with simulation energy as well as the interaction products, it 
is as good as at least within 15% 43|, |45|. 



V. CONCLUSION 

In addition to Galactic 7-ray emission, the Galactic neutrino emission could provide 
complementary information on the existence of hadronic processes in astrophysical objects 
and also on the origin of Galactic cosmic-rays. Here we have calculated the e^, z/g and 
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z/^ production in cosmic-ray interactions and derived matrices that allow to estimate the 
spectra of secondary leptons for arbitrary spectra of the parent hadronic cosmic rays. 

With the production matrices, we calculate the synchrotron radiation of secondary elec- 
trons in the shell-type SNR RX J1713. 7-3649 under the assumption that the observed TeV- 
scale gamma-ray emission is hadronic in origin. We find that the radio synchrotron flux 
of the secondary electrons exceeds the observed flux level if the magnetic field strength is 
too high. The multi-mG fields recently invoked to explain the X-ray flux variations 18| are 
therefore unlikely to extend over a large fraction of the radio-emitting region or the spectrum 
of hadronic cosmic rays in the 0.1-100 GeV energy window must be unusually hard. 

Our calculations for neutrinos show that the neutrino detection rates from the TeV 7-ray 
sources RX J1713. 7-3946 and Vela Jr. are promising for muon energies of a few TeV, or in 
other words, about 10 TeV neutrino energy. An accumulation of data over about five to ten 
years would allow testing the hadronic origin of TeV 7-rays. 

The production matrices for e~, e"*", z/g, Ve-, t'/x ^"^^ ^tJ. ^^^^ t)e made available for download 
at 



http://cherenkov. physics, iastate. edu/lepton-prod 
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FIG. 1: The synchrotron radiation contributed by secondary electrons produced in RX J1713.7- 
3946 assuming a source distance of 1 kpc and the SNR age as 1,600 years. The cosmic-ray particle 
spectrum follows ([8]) and ([9]) with a = 1.8. The curves show the synchrotron spectra for magnetic 
fields B = 500, 2000, and 6000jiG. The ATCA data (open marks) are based the brightness of the 
NW rim of RX J1713. 7-3946 135,1 scaled to the entire remnant with a factor 3 based on the angular 

n n 

distribution of TeV gamma-ray emission [13]. The data points (closed marks) are taken from |13l |. 
where they are referenced as private communication. See text for discussions. 
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FIG. 2: The 7-ray and raw z/^ spectra from the cosmic-ray source RX J1713. 7-3946 derived us- 
ing the production matrix discussed in this work, shown in comparison with results of pubHshed 
parametrizations 
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251 ] based on the same cosmic-ray particle spectrum. With the same 
generating cosmic-ray particle spectrum, our production matrix gives about 15-20% more 7-rays 
at all energies. Our calculation also shows about 30% more neutrinos than the parametric model 
of Kappes et al. 



251] after considering all relevant decay processes in the cosmic-ray hadronic inter- 
actions. The calculations on i/^ spectra indicate that the parametric model of Vissani [22] gives an 
overall lower raw v^ yield for the same 7-ray flux before taking neutrino oscillations into account. 
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FIG. 3: Integrated fully mixed {u^ + t'^)-induced niuon rates from the TeV 7-ray sources 
RX J1713. 7-3946 (solid line), Vela Jr. (dash-dot hne) and MGRO J2019-h37 (dash line) above a 
given muon energy within an IceCube-like detector. The gen erating cosmic-ray spectra are obtained 
jy fitting the observed TeV 7-rays of RX J1713.7-3946 



1$, Vela Jr. 



331, and MGRO J2019-h37 



39[, using the production matrices for gamma-rays and neutrinos. The atmospheric background 
rates, vertical (long dot) and horizontal (short dot), are calculated using the parametrization of 
42l | and assuming an angular resolution element of 3 deg'^. 
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TABLE I: Relative contribution of synchrotron radiation by secondary electrons in RX J1713.7- 
3946 to the observed radio flux in different source models and magnetic fields near the SNR shell. 
The relative contribution of synchrotron radiation in this table are calculated relative to the radio 
data of integrated flux density 20 Jy at 1.4 GHz, which is the result of [35[| scaled to the entire 



remnant with a factor 3 based on the angular distribution of TeV gamma-ray emission 



Q. 



index 


Age (Yr) 


B(/iG) 


ratio 


B (/iG) ratio 


B (/iG) ratio 


B (/iG) ratio 


a = 1.6 


1,600 


500 


3.4% 


1000 8% 


2000 17% 


4000 40% 


a = 1.8 


1,600 


500 


8% 


1000 18% 


2000 41% 


4000 95% 


a = 1.8 


10,000 


500 


48% 


1000 111% 


2000 255% 


4000 594% 


a = 2.0 


1,600 


500 


21% 


1000 50% 


2000 119% 


4000 279% 
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